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ABSTRACT 

The evolution of the cluster mass function and the cluster correlation func- 
tion from z = to z ~ 3 are determined using ~10 6 clusters obtained from 
high-resolution simulations of the current best-fit ACDM cosmology (Q m = 0.27, 
erg = 0.84, h = 0.7). The results provide predictions for comparisons with fu- 
ture observations of high redshift clusters. A comparison of the predicted mass 
function of low redshift clusters with observations from early Sloan Digital Sky 
Survey data, and the predicted abundance of massive distant clusters with obser- 
vational results, favor a slightly larger amplitude of mass fluctuations (<j 8 ~ 0.9) 
and lower density parameter (Q m ~ 0.2); these values are consistent within l-a 
with the current observational and model uncertainties. The cluster correlation 
function strength increases with redshift for a given mass limit; the clusters were 
more strongly correlated in the past, due to their increasing bias with redshift— 
the bias reaches b ~ 100 at z=2 for M > 5 x 10 13 /i -1 M Q clusters. The richness- 
dependent cluster correlation function, represented by the correlation scale versus 
cluster mean separation relation, Ro — d, is generally consistent with observations. 
This relation can be approximated as Rq = 1.7d°' 6 /i _1 Mpc for d ~ 20 — 60 h^ 1 
Mpc. The R — d relation exhibits surprisingly little evolution with redshift for 
z < 2; this can provide a new test of the current ACDM model when compared 
with future observations of high redshift clusters. 

Subject headings: cosmology:observations-cosmology:theory-cosmological parameters- 
dark matter-galaxies:clusters:general-large-scale structure of the universe 

1. Introduction 

Clusters of galaxies, the most massive virialized structures in the universe, provide vital 
information about large-scale structure of the universe and place powerful constraints on 
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cosmology (Bahcall 1988; Peebles 1993; Carlberg et al. 1997; Rosati, Borgani & Norman 
2002; Henry 2004; and references therein). The abundance of clusters as a function of 
mass (i.e. the cluster mass function) and the evolution of this abundance with redshift are 
sensitive probes of both the present day density parameter (fl m ) and the amplitude of mass 
fluctuations (er 8 ); this provides a powerful test of the cosmological model (Peebles, Daly, & 
Juszkiewicz 1989; Henry & Arnaud 1991; Bahcall & Cen 1992; Oukbir & Blanchard 1992; 
Barlett & Silk 1993; White et al. 1993; Viana k Liddle 1996; Eke, Cole, & Frenk 1996; Pen 
1996; Cen 1998; Ikebe et al. 2002; Seljak 2002; Bahcall et al. 2003b; and references therein). 

The spatial distribution of clusters of galaxies serves as a complementary test of the 
cosmological model; this cluster distribution is often described by the two-point cluster 
correlation function. The amplitude of the correlation function offers a strong test of the 
cosmology (Bardeen, Bond, & Efstathiou 1987; Bahcall & Cen 1992; Mann, Heavens & 
Peacock 1993; Holtzman & Primack 1993; Croft & Efstathiou 1994; Borgani et al. 1995). In 
fact, it was the unexpectedly strong observed cluster correlations (Bahcall & Soneira 1983) 
that provided the first evidence against the then standard Q m = 1 SCDM models (Bahcall 
& Cen 1992; Croft et al. 1997; Borgani, Plionis, & Kolokotronis 1999; Governato et al. 
1999; Colberg et al. 2000; and references therein). In addition, the evolution of the cluster 
correlation function with redshift, which has received comparatively less attention in the 
literature, is sensitive to the cosmology. 

Taken together, the cluster mass and correlation functions provide two of the most 
powerful constraints on cosmological models. Predictions have become increasingly robust 
with larger and higher resolution cosmological simulations, made possible by recent growth in 
computing power combined with more sophisticated algorithms. Considerable progress has 
also been made on the observational front to determine cosmological parameters. Combining 
the recent WMAP data with finer-scale CMB experiments plus galaxy and Ly-a forest data, 
Spergel et al. (2003) determined a best-fit power law ACDM model. While further data 
will refine this model, the differences are likely to be small. As a further check of this 
model, we present in this work the simulated mass and correlation functions of clusters of 
galaxies, and their evolution with redshift, determined from mock sky survey cluster catalogs 
generated from a ACDM simulation of the current best-fit cosmological model (Spergel et al. 
2003). We compare our results to the most recent observations, and lay the groundwork for 
comparison with future cluster observations at both low and high redshift. Such comparisons 
will provide important tests of the current cosmology, and will enable further improvements 
in the determination of cosmological parameters. 
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2. The Cluster Mass-Function and Its Evolution 

For the simulation parameters we took those determined by Spergel et al. (2003) from 
the new WMAP data on the largest scales, supplemented by other CMB experiments, galaxy 
surveys, and Ly-a forest observations on smaller scales. Assuming a spatially flat power law 
ACDM universe, these are: matter density Vl m = 0.27, cosmological constant A = 0.73, 
power spectrum amplitude cr 8 = 0.84 , spectral index n s = 0.96, and h = 0.7 where H = 
lOO/i km-s^-Mpc -1 . The simulation used the TPM code (Bode & Ostriker 2003) to evolve 
N = 1260 3 « 2 x 10 9 particles in a periodic box 1500/i 1 Mpc on a side. The particle mass 
is m, p = 1.26 x lO n /i -1 M , and a cubic spline softening length of 17 fo _1 kpc was introduced. 
The simulation is discussed in more detail in Hopkins, Bahcall & Bode (2005). 

Particle positions in a light cone covering one octant of the sky out to redshift z = 3 
were saved to disk; snapshots of the entire simulation volume were also saved. Dark matter 
halos, which would house clusters of galaxies, were identified from the particle position data 
using a Friends-of-Friends (FOF) percolation algorithm, with a linking length of b = 0.2 
times the mean particle separation (Lacey & Cole 1994). The cluster center was defined as 
the location of the most bound particle. Clusters identified using linking length parameters 
of b = 0.16 and 0.25 were examined for comparison, yielding similar results. 

The mass function (MF) of clusters, n(>M), represents the number density of clusters 
with mass greater than M. The constraints which the present day cluster MF places on the 
mean density parameter of the universe (fl m ) and the amplitude of mass fluctuations (cr 8 ) are 
partially degenerate in f2 m — er 8 . Observations of the present day cluster MF have established 
that a 8 tt° m 5 « 0.5 (Henry & Arnaud 1991; Bahcall & Cen 1992; White, Efstathiou, & Frenk 
1993; Eke, Cole, & Frenk 1996; Viana k Liddle 1996; Kitayama & Suto 1997; Pen 1998). This 
degeneracy can be broken by studying the evolution of the cluster MF with redshift (Peebles, 
Daly, & Juszkiewicz 1989; Oukbir & Blanchard 1992, 1997; Eke, Cole, & Frenk 1996; Viana 
& Liddle 1996; Bahcall, Fan, & Cen 1997; Carlberg, Morris, Yee, & Ellingson 1997; Henry 
1997, 2000; Bahcall & Fan 1998; Eke et al. 1998; Donahue & Voit 1999). Cluster evolution 
is exponentially dependent on cr| (Bahcall, Fan, & Cen 1997; Bahcall & Bode 2003), and as 
a result, the most massive clusters evolve strongly in a low-erg, VL m = 1 universe, producing 
a very low abundance of massive clusters at z > 0.5. Conversely, the evolution of rich 
clusters is significantly weaker in a cr 8 1 \ow-Q m universe, with a considerably higher 
cluster abundance at z > 0.5 as compared to lower-a 8 models. 

The simulated cluster MF was determined from the light cone outputs using cluster 
masses calculated according to typical masses used by observers, including: mass within 
fixed radii (relative to the center of the cluster) of 0.5 h^Mpc comoving (M 5 ), 1.5 h^Mpc 
comoving (M L5 ), 0.6 h^Mpc physical (M a6 ), and also mass within a radius containing a 



-4- 



mean overdensity of 200 relative to the critical density (M 2 oo)- Minimum mass cutoffs were 
chosen in order to ensure the completeness of the cluster sample: M 05 > 1.6 x 10 13 /i _1 M Q ; 
Mi.5 > 5 x 1O 13 /i _1 M ; M . 6 > 3 x 1O 13 /i _1 M ; and > f .75 x 10 13 /i _1 M o . The evolution 
of the cluster MF for M .5 is presented in Figure 1. These results can be used for comparison 
of predictions of the current cosmological model with future observations of high redshift 
clusters. 

The cluster MF for M 6 and z — 0.1 — 0.2 was compared to the observed early Sloan 
Digital Sky Survey cluster MF for the same redshift and mass range (Bahcall et al 2003b; 
see also Bahcall et al. 2003a). The results are presented in Figure 2; also shown are the 
best analytic model fits. The observed data follow the same shape as the LCDM MF, but 
are systematically offset to slightly lower masses. This suggests that either a bias of ~20% 
(~l-sigma) exists in the observed cluster mass calibration (the mass calibration is expected 
to improve for the larger upcoming sample of SDSS clusters) or that Q m is somewhat lower 
than used in the simulation. The best-fit parameters to the observed SDSS MF are £l m ~ 0.2 
and as ~ 0.9, as shown in Figure 2. These best-fit parameters are consistent with the current 
model parameters within the combined observational and model uncertainties. 

The predicted abundance evolution of high mass clusters can be compared to observa- 
tions at higher redshift. We use the mass threshold Mi. 5 > 8 x 1O 14 /i _1 M and observed 
abundances of Bahcall & Bode (2003); the results are presented in Figure 3. At high red- 
shift, the model predicts considerably fewer clusters than observed. This suggests that the 
amplitude of mass fluctuations a s is larger than 0.84. Bahcall & Bode (2003) found a best 
fit value of as = 0.9 ±0.1 for these data. These results are consistent with the current 
cosmological parameter (cr 8 = 0.84 ±0.04) within the combined observational and model un- 
certainties. Both the local cluster MF and the evolution of massive clusters at high redshift 
thus suggest <t 8 w 0.9 and tt m ps 0.2. These results indicate a a 8 value at the high end of 
the current best-fit model allowed range, with Vt m at the lower end of the accepted range; 
both values are consistent within 1-sigma with the recent CMB and large-scale structure 
observations. Further observations are needed to narrow down the precise value of these 
cosmological parameters. 



3. The Cluster Correlation Function and Its Evolution 

The cluster correlation function (CF) is a statistical measure of how strongly clusters of 
galaxies cluster as a function of scale. The probability of finding a pair of clusters in volumes 
V\ and V2, as a function of pair separation (r) is 

dP = n 2 {l + i cc {r))dV 1 dV 2) (1) 
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where n is the mean number density of clusters, and £ cc (r) is the cluster CF. The spatial 
distribution of clusters of galaxies described by the cluster CF is sensitive to cosmological 
parameters (e.g. Bahcall & Cen 1992; Borgani, Plionis, & Kolokotronis 1999; Colberg et al. 
2000; and references therein). 

Observationally, the cluster CF is an order of magnitude stronger than that of individual 
galaxies: typical galaxy correlation scales are ~ 5/i _1 Mpc, as compared to ~20 — 25/i _1 Mpc 
for the richest clusters (Bahcall & Soneira 1983; Klypin & Kopylov 1983; see also Bahcall 
1988; Huchra et al. 1990; Postman, Huchra, & Guller 1992; Bahcall & West 1992; Peacock 
& West 1992; Dalton et al. 1994; Croft et al. 1997; Abadi, Lambas, & Muriel 1998; Lee & 
Park 1999; Borgani, Plionis, & Kolokotronis 1999; Collins et al. 2000; Gonzalez, Zaritsky, & 
Wechsler 2002; and references therein). Furthermore, the strength of the CF increases with 
cluster richness and mass (Bahcall & Soneira 1983). As a result, the rarest, most massive 
clusters exhibit the strongest correlations. 

The richness dependence of the cluster CF has been confirmed observationally (see 
references above), and explained theoretically (Kaiser 1984; Bahcall & Cen 1992; Mo & 
White 1996; Governato et al. 1999; Colberg et al. 2000; Moscardini et al. 2000; Sheth, 
Mo, & Tormen 2001; and references therein). However, these analyses have been done 
at low redshift, z < 0.5. With observational data becoming available at higher redshifts, 
the expected evolution of the cluster correlation function is increasingly important as an 
independent test of the cosmological model. Analytic approximations to the evolution of 
cluster halo abundance, bias, and clustering have yielded some promising results (Mann, 
Heavens, & Peacock 1993; Mo & White 1996, 2002; Sheth, Mo, & Tormen 2001; Moscardini 
et al. 2001). However, numerical simulations can provide the most reliable comparison 
between theory and observations. We determine the expected cluster CF and its evolution 
for the current best-fit ACDM model using light-cone outputs from the N-body cosmological 
simulation discussed in §2. When compared with recent observational results at low redshift, 
this provides a test of the current cosmological model; at the same time the evolution of the 
ACDM CF provides detailed predictions for comparison with future observations of high 
redshift clusters. 

The cluster CF was calculated as a function of separation using £ cc (r) = F DD {r)/ F RR {r) — 
1, where Fdd{t) is the frequency of cluster pairs with comoving separation r, and F RR {r) is 
the frequency of pairs in a random catalog. The cluster CF was calculated for different mass 
thresholds for M . 5 and M 200 in several redshift bins: z= 0-0.2, 0.45-0.55, 0.9-1.1, 1.4-1.6, 
and 1.9-2.2. Because of the rapidly decreasing abundance of the most massive clusters with 
redshift, their CF is studied only at lower redshifts. Examples of the cluster CF at different 
redshifts are presented in Figure 4. 
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The cluster CF for each redshift bin and mass threshold was fit to a power law of the 
form £ cc (r) = {r/R )~' y , where 7 is the correlation slope and Rq is the correlation scale. The 
fits were done over the linear range of the cluster CF (r < 50/i -1 Mpc) for both a fixed slope 
of 7 = 2, and for 7 as a free fitting parameter. The results are similar for both methods, 
with the best fit free slope 7 ~ 1.8 for z < 1.5. The slope is mildly richness-dependent, 
with more massive clusters showing a slightly steeper slope. The evolution of Rq and 7 with 
redshift is presented in Figure 5. 

The correlation scale Rq increases both with cluster mass and with redshift (see Figure 
5). The steepening slope at high redshifts causes Rq{z) to be slightly lower for a free 7 
as compared to that of a fixed 7 = 2, but the evolutionary trend remains the same. The 
evolutionary increase of the cluster correlation scale with redshift is stronger for more massive 
clusters at higher redshift. For example, clusters with M .5 > 1.6 x lO 13 /i _1 M have Rq 
10/i _1 Mpc at z = 0.15, R « 13/i -1 Mpc at z = 1, and Rq « 17/i _1 Mpc at z = 2; while 
clusters with M . 5 > 3.0 x 10 13 /i" 1 M o have R « 12/i -1 Mpc at z = 0.15, R « 17/i _1 Mpc at 
z = 1, and i? ~ 29/?r 1 Mpc ai z — 2. The free slope fits yield similar results. 

The cluster CF (using M FO f) and the CF of the dark matter particles were also deter- 
mined for simulation box snapshots at various redshifts, and fit to a power law with a fixed 
slope 7 = 2. The correlation scale of clusters {Rq), as before, increases with redshift. By 
contrast, the correlation scale of the mass {Rq 1 ) decreases with redshift, as expected. This 
is due to the fact that clusters of a given comoving mass represent higher density peaks of 
the mass distribution as the redshift increases, and thus exhibit enhanced clustering with 
increasing redshift (see also Cole & Kaiser 1989; Mo & White 1996, 2002; Sheth, Mo, & 
Tormen 2001; Moscardini et al. 2001). Figure 6 shows the evolution of the ratio Rq/Rq 1 
with redshift, which follows the evolution of bias from z = to 2, where bias is defined as 
b 2 = icc/imm ~ {Rq/Rq 1 ) 1 - As seen in Figure 6, the ratio Rq/R™, and thus the bias b, 
increases from ~3 to ~100 (for M FOF > 5 x lO 13 /i _1 M clusters) as the redshift increases 
from z = to 2. Other mass clusters show a similar trend. 

Another useful approach to studying the evolution of the cluster correlation function is 
the Rq — d relation, where Rq is the fitted correlation scale and d is the mean intercluster 
comoving separation. For a larger mass limit objects are less abundant, and thus their mean 
separation d increases. These objects are also more biased, so an increasing Rq with d is 
observed (Bahcall & Soneira 1983; Szalay & Schramm 1985; Bahcall 1988; Croft et al. 1997; 
Governato et al. 1999; Bahcall et al. 2003c; Padilla et al. 2004; and references therein). 
Thus, the evolution of the Rq — d relation allows us to investigate the change in correlation 
strength with cluster mass and with redshift. We present the Rq — d relation (using M , 5 
mass thresholds and a fixed 7) for several redshift bins in Figure 7. The results for a free 7 
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are shown in Figure 8. 

The resulting R —d relation shows a surprising behavior; there is essentially no evolution 
with redshifts for z < 2 (using M . 5 thresholds). Mass thresholds measured within a fixed 
overdensity (M 20 o and Mfof) show slightly more evolution, but still surprisingly little. This 
redshift invariant Rq — d relation provides a powerful test of the cosmology when compared 
with upcoming observations of high redshift clusters. 

Why is the Rq — d relation invariant with redshift? The invariance appears to be partly 
due to the relative constancy of the cluster mass hierarchy with redshift. That is, the halos 
which are the most massive at an early time tend to remain part of the population of the 
most massive objects. Consider a given comoving volume of space. The majority, but not all, 
of the N most massive clusters at z = are among the N most massive clusters at higher 
redshift. The clusters that are not in the top N at higher redshift are, in turn, clusters 
lower down on the mass scale. In such a case, if the sample of clusters at a given mean 
comoving separation d is not dramatically different at different redshifts, it will yield similar 
correlation strengths R , assuming that the clusters have not moved significantly over that 
time period. To confirm this, we select clusters using M 05 from box snapshots at z—0, 
0.94, 1.4, and 2.0. A cluster at high redshift is considered to be the 'matched' system if it 
is within 3 h~ 1 Mpc physical separation of its position at z — 0. Using a fixed comoving 
mean separation d = (N/V) -1 / 3 at all redshifts, the N m clusters that match with their z = 
counterpart are kept, and the remaining N — N m clusters are selected randomly from the next 
N clusters down the mass ladder. Using this new distribution, we find that the correlation 
scale Ro(d, z) is indeed nearly constant with z for a given d (for d < 90/i _1 Mpc and z < 2). 
When looking along the past light cone, the Rq — d relation is thus nearly independent of 
redshift for z < 2 because the majority of clusters at different redshifts represent a similar 
population along the same filamentary hierarchical structure, near the top of the cluster 
mass hierarchy. 

We compare the predictions of the ACDM model with current Rq — d observations in 
Figure 9, using M a5 thresholds and a fixed correlation slope of 7 = 2 (for details on the 
observations see Table 1 of Bahcall et al. 2003c). All the observed correlation scales (R ) 
and mean separations (d) were converted to comoving scales in a ACDM cosmology. The 
band in Figure 9 represents the simulated Rq — d relation (with 1-cr range) at z = — 0.3 
(to match the redshift range of the observations). An analytic approximation to the ACDM 
Rq — d relation for 20 < d < 60/i _1 Mpc, presented by the dashed curve, is 



The results show a general agreement between the ACDM model and observations. The op- 
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tically selected cluster samples agree with the model within 1-er. There is, however, a wide 
scatter at the high-ci end, especially when X-ray selected clusters are included; the X-ray 
selected clusters seem to suggest a somewhat larger R than optically selected clusters. Due 
to this large scatter, the correlation function cannot yet provide a high-precision determina- 
tion of cosmological parameters. It does, however, clearly rule out high-f2 m models, as their 
far weaker correlations are inconsistent with the data. The strong correlations suggested by 
the current X-ray clusters will be tested with future observations of X-ray selected cluster 
samples; this should clarify whether or not the X-ray clusters are consistent with the optical 
data and with the best-fit ACDM cosmology. 



4. Conclusions 

We use light cone outputs from a large-scale simulation of the currently best-fit ACDM 
cosmological model to generate mock sky survey cluster catalogs that can be readily com- 
pared with observations. The catalogs were used to determine the present day cluster mass 
and correlation functions, which together constitute a sensitive test of the cosmological 
model. We determine the evolution of the cluster mass function from z = to z = 3, and 
the evolution of the cluster correlation function from z = to z = 2. These results provide 
predictions of the current cosmological model for comparison with future observations of 
high redshift clusters. Such comparisons can be used to test the current model and provide 
new and independent constraints on both the cosmological density parameter, Q m , and the 
amplitude of mass fluctuations, a$. 

The simulated cluster mass function at low redshift (z—0. 1-0.2) is compared with the 
early Sloan Digital Sky Survey cluster mass function (SDSS: Bahcall et al. 2003b). The 
ACDM mass function predicts somewhat higher abundances than are observed. This sug- 
gests that either a small bias (~20%) exists in the observed cluster mass calibration, or a 
somewhat lower value of the cosmological density parameter is needed (Q m « 0.2). The re- 
sults, however, are consistent with the current cosmology within the combined observational 
uncertainties. The evolution of the most massive clusters is exponentially dependent on cr| 
(Bahcall & Fan 1998), and therefore can be used to break the Q m — os degeneracy that 
exists in the low-redshift mass function. We find that the ACDM model predicts a consid- 
erably lower abundance of distant (z > 0.5) massive clusters than observed. This suggests a 
normalization of <7$ ~ 0.9. 

We determine the cluster correlation function in the ACDM cosmology as a function 
of cluster mass threshold, summarized in the Rq — d relation (Figure 7-8). We find good 
agreement with observations of optically selected clusters (Figure 9); X-ray selected clusters 
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appear to suggest somewhat increased cluster correlations. The wide scatter in the obser- 
vational data at high-rf, makes it difficult at this time to provide precise constraints on the 
cosmology. The data are, however, in good agreement with the current ACDM model. 

We determine the evolution of the richness dependent cluster correlation function in the 
ACDM cosmology (Figures 4 — 8). For a given mass limit, the correlation strength increases 
with redshift. Surprisingly, the Rq — d relation shows no significant evolution out to z < 2. 
This surprising result provides a new, independent test of the current cosmological model 
when compared with future observations of high redshift clusters. 

We thank J. P. Ostriker, D. Spergel, and M. Strauss for helpful discussions. The com- 
putations were performed on the National Science Foundation Terascale Computing System 
at the Pittsburgh Supercomputing Center. This work was supported in part by NSF grant 
AST-0407305. 
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5xl0 13 10 14 
M(<0.5 h" 1 Mpc Comoving) h" 1 (M ) 

Fig. I. — The simulated cluster mass function for ACDM, in redshift bins z— 0.0-0.15; 0.45- 
0.55; 0.95-1.05; 1.45-1.55; 1.9-2.1; 2.2-2.5; and 2.5-3.0 (top to bottom). The mass M . 5 is 
measured within a fixed comoving radius 0.5 ft, _1 Mpc from the cluster center. 
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Fig. 2. — A comparison of the SDSS observed cluster MF and analytic best fits for the HMF 
(filled square, dashed line) and maxBCG (open square, dotted line) selection techniques 
(Bahcall et at. 2003a,b) and the current ACDM simulation predictions for z = 0.1 —0.2. All 
masses are defined within a physical radius of 0.6 h^ 1 Mpc of the cluster center. 
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Fig. 3. — The cluster abundance evolution for massive clusters (M1.5 > 8 x 10 14 h 1 M ) in 
the ACDM model (dashed line) versus observed abundances (1-er Poisson error bars). 
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Fig. 4. — The simulated cluster correlation function for different redshift bins From bottom 
to top: filled circles, z=0.0-0.2; open circles, 0.9-1.1; filled triangles, 1.4-1.6; open squares, 
1.9-2.2; with 1-er Poisson error bars. The separation (r) is in comoving units, and the masses 
are measured within a fixed comoving radius of 0.5 h^ 1 Mpc. The lines represent the best 
fit power laws for a fixed correlation slope 7 = 2. 
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Fig. 5. — Evolution of the cluster correlation scale (Rq comoving) and slope (7) as a function 
of redshift for different mass thresholds for the ACDM model. 
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Fig. 6. — The ratio of the correlation scale of clusters (Rq) to that of the underlying mass 
distribution (Rq 1 ) as a function of redshift for different mass thresholds. Also shown are the 
evolution of R™ (dashed line) and Rq, separately. 
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Fig. 7. — The evolution of the richness-dependent cluster correlation function : the Rq — d 
relation (comoving coordinates) for different redshifts using a fixed correlation slope 7 = 2. 
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Fig. 8. — The R — d relation (comoving coordinates) for different redshifts using a free- fitting 
correlation slope (7). 
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Fig. 9. — A comparison between the ACDM model Rq — d relation (shaded band; z=0-0.3) 
and observational results for a fixed correlation slope 7 = 2 (data compiled by Bahcall et al. 
2003c, with the addition of the 2PIGG groups from Padilla et al. 2004). The dashed line is 
a best fit power law approximation: R = l.ld°- & h~ l Mpc. 



